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Abstract 

We consider different generalizations of the Fokker-Planck-equation devised to de- 
scribe Levy processes in potential force fields. We show that such generalizations 
can proceed along different lines. On one hand, Levy statistics can emerge from the 
fractal temporal nature of the underlying process, i.e. a high variability in the rate 
of microscopic events. On the other hand, they may be a direct consequence of the 
scale-free spatial structure on which the process evolves. Although both forms con- 
sidered lead to Boltzmann equilibrium, the relaxation patterns are quite different. 
As an example, generalized diffusion in a double-well potential is considered. 



1 Introduction 



Random walk processes leading to anomalous diffusion are adequate for de- 
scribing various physical situations. The continuous time random walk (CTRW) 
model of Scher and Montroll [1] for instance, leading to strongly sub-diffusive 
behavior, was a milestone in understanding photoconductivity in strongly dis- 
ordered and glassy semiconductors. Due to its simplicity, this model was re- 
cently employed to investigate aging phenomena typical for glasses and re- 
lated complex systems [2,3,4,5]. On the other hand, Levy-flight models [6], 
leading to superdiffusion, are adequate for the description of phenomena rang- 
ing from transport in heterogeneous catalysis [7], self-diffusion in micelle sys- 
tems [8], reactions and transport in polymer systems under conformational 
motion [9], transport processes in heterogeneous rocks [10] and the behavior 
of dynamical systems [11] to flight paths of albatrosses [12] and even human 
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eye-movements [13]. Closely related models appear in the description of eco- 
nomic time series [14]. Levy-related statistics were observed in hydrodynamic 
transport [15], and in the motion of gold nanoclusters on graphite [16]. In addi- 
tion, mixed models were proposed in which slow temporal evolution (described 
by Scher-Montroll CTRW) is combined with the possibility of Levy-jumps, so 
that in general both sub- or superdiffusive behavior can arise [17,18]. Similar- 
ily, processes often referred to as Levy- walks are based on long-tailed jumps in 
combinations with a time cost to perform them, see for example [6,19]. Levy- 
walks were introduced in [20] to explain accelerated diffusion in Josephson 
junctions and related chaotic systems [21]. 

Typically, the corresponding processes are described on a stochastic level, such 
that the incorporation of boundary conditions, time- or position-dependent 
forces leads to considerable mathematical difficulties. 

Concerning ordinary diffusion and transport in a stochastic system, the Fokker- 
Planck equation (FPE) is an adequate instrument for solving boundary value 
or time-dependent problems [22]. The FPE is an equation for the probability 
density function (pdf ) p(x, t) of the particle's position at time t and reads 

= C FPP = -V (fif(x,t)p(x,t)) + K Ap(x,t). (1) 

Here Cpp denotes the Fokker-Planck operator, \i and k are the mobility and 
the diffusion coefficient (assumed to be constant) and f(x,t) is the external 
force acting on the particle. Note that since the mobility and the diffusion 
coefficient are connected via the Einstein relation, k = kTji = fx//3, the equa- 
tion effectively contains only one free parameter. Alternatively, eq. (1) may be 
recast into a different form in the case of a potential force f(x, t) = — V$(x, t): 



dp(x, t) 
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(2) 



the equivalence of which to eq. (1) is easily checked by expanding the differ- 
entiations. Obviously, the stationary solution of the FPE (if it exists) is the 
Boltzmann distribution, p(x) oc exp(— /3<E>). Inserting this form into eq. (2) one 
readily infers that the two terms entering with opposite signs are equal and 
cancel. 



In order to investigate analogous systems incorporating anomalous diffusion, 
various generalizations of the FPE were proposed. Concerning subdiffusion 
the appropriate generalizations can be obtained from CTRW models, for a 
general survey see [27], as well as from an approach based on fractional master 
equations [28]. Concerning superdiffusion (Levy-flights), such generalizations 
typically involve changing from a gradient and Laplacian to the corresponding 
fractional derivative. For example, a pure Levy-flight can be described by 
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means of the fractional diffusion equation [23,24,25] 

^M = « aA ^p(x,0, (3) 

where A"/ 2 is a fractional Laplacian, a linear operator, whose action on a 
function f(x) in Fourier space is described by A a ^ 2 f(x) = — (k 2 ) a ^ 2 f(k) = 
— \k\ a f(k). The coordinate representation of this operator is discussed in Ap- 
pendix 1. The solution of eq. (3) in Fourier space reads: 

G(k,t) =exp{-K\k\ a t) (4) 

which in coordinate space corresponds to a symmetric Levy stable distribution 

G(x ' t) = (^ L ((^ ;a '°) (5) 

(we use here the canonical notation, see ref. [26]). On the other hand, there 
is no generic way to generalize the Fokker-Planck operator containing an ex- 
ternal force term: For example, the drift term may stay unchanged, may cor- 
respond to a symmetric or asymmetric fractional derivative, or the whole 
Fokker-Planck-operator may be raised to the power of a/2. All correspond- 
ing equations have their physical meaning and applications. The properties 
of their solutions differ considerably. Moreover, as we proceed to show, other 
generalizations are also possible. The situation here is to some extent simi- 
lar to one in fractal geometry: What is the generalization of the Euclidean 
dimension to the fractal case? Depending on the system and on the prop- 
erty of interest it could be the fractal dimension, the spectral dimension, or 
even a spectrum of different dimensions. In what follows we discuss the prop- 
erties of some generalization of the FPE to the superdiffusive case leading 
to Levy-flights in the force- free limit. We will focus on two complementary 
situations, namely one in which superdiffusion stems from fractal temporal 
properties of the underlying microscopic dynamics and another in which the 
spatial structure is responsible for enhanced diffusion. We mostly concentrate 
on the generalizations of the forms, eq. (1) and eq. (2) which lead to a Boltz- 
mann distribution in equilibrium. We discuss the physical implications of the 
corresponding equations and some properties of their solutions. As a physically 
most interesting application, a fractional generalization of a Kramers problem 
in a double-well potential is considered. Let us first turn to the situations 
which can be considered as stemming from the fractal temporal behavior. 



2 Levy processes as stemming from subordination 

An important property of symmetric Levy processes is the fact that they 
are subordinated to ordinary Brownian motion: the corresponding pdf can be 
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represented in the form: 

P(*. t) = ±- a L ( JL_; a, 0) = / -±= exp ( JL; ° l) dr 

(6) 

where L (x; a/2, 1) is an extreme (one-sided) Levy distribution of index a/2 [26]. 

The variable r is called the operational time of the process. The interpretation 
of eq. (6) is that Levy-flights can be considered as stemming from a highly 
irregular sampling trajectories generated by simple diffusion (random walk). 
This type of temporal behavior plays a very important role in many situations 
encountered in physics, such as anequilibrium phase transitions [29]. For a 
review see [30] and references therein. 

Here, the trajectory of diffusion (a random walk in a discrete case or a Wiener 
path in a continuous situation) is parameterized by the operational time r (say, 
the number of steps of the random walk), which itself is a random function 
of the (physical) time. The random process r(t) is a process with positive 
increments, and the distribution of r(t) is given by a one-sided Levy law. 
Thus, the increase At of the operational time per physical time unit is given 
by a one-sided Levy-distribution, p(Ar) = L (At; a/ 2,1) having a power-law 
tail, p(Ar) oc t- 1 -^ 2 . 

The corresponding situation is depicted in fig. 1. Here a Levy-flight of N — 100 
steps is generated as follows. One generates a sequence {rii} (i = 1,2,..., N) of 
the numbers of steps between the observation instants and plots a realization 
of a simple random walk of M — Y.f=i n i steps. The turning points of a Levy- 
flight are then the positions of a random walker after rrik steps, = Yh=i n i- 
Fig. 1 corresponds to a process subordinated to a random walk under the 
operational time given by 

p( " ,) = ^F exp (-i)- (7) 

The overall pdf of the process p(x, t) is then given approximately by a Cauchy 
distribution, 

p(a; ' t)= S7s exp (-97^ exp ("£) " 

"/7t exp ("97^ exp ("£)* I = ^T^)' (8) 

Although this continuous approximation is valid here, some care must be taken 
in the general case when sums over steps in random walks are replaced by 
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Figure 1. The mechanism of subordination: Levy-flights can be considered as stem- 
ming from a (highly inhomogeneous) random sampling of a simple random walk. 

integral, see for example [31]. Applying the subordination procedure to the 
general, biased diffusive process, one is lead to the generalization of the FPE, 
in which the entire Fokker-Planck-operator is raised to the power a/2: 

— = - {-C F p) p (9) 

as shown in [32]. The generalized Fokker-Planck operator C a = — {—Cpp) 01 ^ 2 
commutes with Cpp and shares with it the same set of the eigenfunctions 
4>i. The corresponding eigenvalues A$ of the operator C a are connected with 
those of Cfp via Aj = — (— Aj) a / 2 . Since the eigenvalues Aj of Cfp are real and 
nonpositive, this is true for the Aj as well. Eq. (9) follows from the subordina- 
tion procedure just in the same way as a Cauchy distribution followed from a 
Gaussian distribution in the example given above. Let us consider a spectral 
representation of the solution of the FPE, 

g(x,r)=Y,Mx)e-^ T . (10) 
The subordination procedure applied to this solution leads to 
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p{x,t) = Jg(x,T)^L^^;^,l ) jdT= (11) 

i 

due to a well-known property of the Laplace-transforms of one-sided Levy dis- 
tributions [26]. Note that the eigenfunction O corresponding to a zero eigen- 
value (if any) is a stationary distribution, which is the same (i.e. a Boltzmann 
distribution) for the normal and for the generalized process. The spectral de- 
composition of solutions of eq. (9) given by eq. (11) is incorporated in our 
numerical procedure for the solution of superdiffusive FPE in a potential field 
considered in section 4. 



Eq. (9) can also be interpreted within the Langevin scheme, which clarifies 
the circumstances under which Boltzmann statistics may or may not appear 
when superdiffusive generalizations of the FPEs are introduced. Let us first 
consider a purely diffusive situation without external force. In the continuous 
limit (corresponding to averaging over the time periods which are short enough 
to consider all parameters constant but during which many steps of a random 
walk are performed), the system's development in its operational time is given 
by a Langevin equation 

d -x(r) = ^(r) (12) 



dr 

where £(r) is a 5-correlated Gaussian noise with zero mean and with unit 
dispersion. Subordination requires that 

Tt = A(t) ' (13 » 

where X(t) is a one-sided Levy process of index a/2. The resulting process is 
thus described by a Langevin equation 

! = ££*M = ^(*««)>. (14) 

where r(t) is a monotonously increasing random function which is chang- 
ing slowly on the scale on which £(r)is correlated. This suggests that the 
random process £(t(£)) is a Gaussian one, has zero mean and dispersion 
(ar(tm(r(f))) = 5 (r(t) - r(f)) = {dr/dt)^5{t - f) = A(t)- 1 ^ - f), 
indicating that the process x(t) can be formally represented by a Langevin 
equation 

X(t) = yj2K\(t) V (t) (15) 

where rj(t) is a 5-correlated Gaussian noise with zero mean and unit dispersion. 
From this equation, a FPE of the type of eq. (3) follows along the usual lines. 
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Physically, such strongly inhomogeneous behavior can be attributed to random 
fluctuations of the mobility or of the temperature or of the both: K,(t) = 
n\(t) = kfi(t)T(t). Let us return to the full process described by a Langevin 
equation 



x(t) =fif + V / WCW (16) 

with a potential force f(x,t) = —VU(x,t). Here the situations corresponding 
to the fluctuations in temperature and to ones in mobility differ vastly. Keeping 
the temperature constant and letting \x fluctuate leads to a process subordi- 
nated to a biased random walk, which is described by a Langevin equation 



x(t) = fi(t)f + ^2fi(t)kT^(t) (17) 

where n(t) is now a strongly fluctuating random function, whose distribution 
is given by a one-sided Levy law. The corresponding pdf is governed exactly 
by a fractional Fokker-Planck equation, eq. (9) derived above. This equation 
shows a relaxation to a normal Boltzmann distribution, which is not surprising 
since the corresponding scheme describes a system showing detailed balance 
at given, fixed temperature T. Note that although the equilibrium properties 
of the system are usual, its relaxation properties are not, see ref. [34] for a 
detailed discussion. Thus, it does not exhibit linear response to a constant 
external field, and, moreover, since the corresponding pdf neither possesses 
the second moment nor scales, no reasonable generalization of the Einstein's 
relation exists. The same is valid for the model of topological superdiffusion 
discussed below. On the other hand, as we proceed to show, the behavior of 
both models in confining potentials such as a harmonic or a double-well one, 
is not as exotic (vide infra). 

The situation under fluctuating temperature is vastly different: a system with 
fluctuating temperature must not show Boltzmann statistics at equilibrium (if 
any) . Let us mention that under certain circumstances the distributions given 
by Langevin equations with fluctuating parameters approach ones given by 
Tsallis statistics [35], i.e. ones stemming from a thermodynamic construct in 
which the entropy is non-extensive and thus the concept of temperature can 
not be introduced in a usual way. 

In our case, fluctuations in temperature do not affect the first term, i.e. the 
deterministic motion and lead (after performing usual steps) to a fractional 
Fokker-Planck equation (FFPE) in the form: 

| = _| Wp)+K A«/ V (18) 

whose solutions are discussed in ref. [36]. 
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Figure 2. Random walk of a particle on a heterogeneous polymer subjected to fast 
conformational changes. The polymer consist of different types of monomers denoted 
by the symbols along the chain (top). Each monomer is associated with an intrinsic 
potential Due to conformational changes, the walker may jump between two 

sites nearby in Euclidean space but far apart along the chemical sequence of the 
polymer. 

3 Topologically induced superdiffusion 

Complementary to the situation discussed in the previous section in which 
superdiffusion originates from fractal temporal behavior, we will discuss the 
fractional generalization of the FPE adequate for the description of diffusion 
on scale-free structures. Contrary to the situation investigated in the previous 
section in which superdiffusion is caused by a fractal sampling of an otherwise 
continuous path, we will focus on the limiting behavior of a master equation 
describing a Markovian jump process in which transitions are caused by ther- 
mal hopping between energetically different states on a topologically complex 
structure, ref. [9]. 

As a model system, we consider a heteropolymer in solution depicted in fig. 2. 
We focus on the dynamics of a test particle performing thermally activated 
motion along the linear chemical coordinate axis x of the polymer. If the 
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polymer is fixed at full extend, thermal activation causes the particle to hop 
at a rate 1/td from a given site to a neighboring site with a probability 
determined by the energy difference within a given mononer pair. However, a 
polymer in solution is subjected to thermal conformational changes occurring 
on a time scale tq. The conformational changes allow sites which are far apart 
along the chemical sequence to come close in Euclidian space, as indicated in 
figure 2: For simplicity, a Gaussian chain (corresponding to a Rouse dynamics) 
can be considered, moreover, the conformational changes will be taken fast 
enough on the scale td- In this case the particle may jump from a given site x 
to a site y far apart with a rate proportional to the probability that the sites 
x and y have come close in the given time interval. In our case, the probability 
of such event is prortional to \x — y\~ 3 / 2 . As we proceed to show, the fractional 
generalization of the Fokker-Planck equation corresponding to the situation 
depicted above differs from the ones considered so far. 

The system introduced above can be modeled by a Markovian jump process 
governed by a temporally homogeneous master equation 



dpjx, t) 
dt 



j iy[w{x \ y)l>M - w( y\z)M (19) 



in which the rates w(x\y) are determined by the short time behavior of the 
conditional probability to make a transition y — > x 

w{x\y) = Urn -^p(x, t + At\y, t). (20) 

Let us mention here, that in eq.(20) the rates w(x\y) may be singular and 
may even defy normalizability. The Cauchy process for example is defined by 
w(x\y) = l/ir(x — y)~ 2 . If, on the other hand, / dxw(x\y) = const, exists the 
constant can be factored out and absorbed by a rescaling of the time. In this 
case the process is a pure jump process, consecutive jumps are seperated by 
time periods during which the random process is constant, in contrast to the 
Cauchy process realizations of which are discontinuous almost everywhere. 

In our example the coordinate x represents the chemical coordinate along the 
polymer chain. The reasoning above suggests a rate of the form 

w (x\ y ) = i e -fl*(*)-*G/)]/2 f( x _ y y (21) 

To 

The right hand side of eq. (21) consists of the thermal component exp(— f3[$(x)- 
$>{%/)]/ 2) accounting for the fact that transitions between energetically differ- 
ent states are less likely to occur if the potential difference <&(x) — $(y) is high. 
The second factor, f(x — y) represents the probability of states with chemical 
coordinates x and y to be nearby in Euclidean space, i.e. this term accounts for 
the geometrical complexity of the system. Due to the translational invariance 
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and symmetry of the system in Euclidean space this term is symmetric and 
depends on the distance \x — y\ only. Obviously, a system defined by a rate 
(21) fulfills detailed balance, characteristic of systems in thermal equilibrium. 
Note also that the stationary solution of Eq. (19) with the rates given by (21) 
is the Boltzmann distribution 

p s (x) oce-^* (x) . (22) 



It is instructive to consider a fully extended polymer of total length L. In this 
case only nearest neighbor hopping occurs, i.e. 

f{x) = l -{5{x-a) + 5{x + a)). (23) 

In (23) cr C L is the spacing between adjacent monomers. Setting qt(x) = 
e l3 ^ x ^ 2 p(x,t) and s(x) = e~^ x ^ 2 the master equation (19) reduces to: 



= ^rs(x) {qt(x + *)+ q t (x - a)} (24) 

—^—q t (x){s(x + a) + s(x-a)} (25) 

_ f v / qt(x + a) + g t (x - a) - 2q t (x) 
~2t S[X) \ 

a 2 | s(x + a) + s(x — a) — 2s(x) 



-s(x) { (26) 



2r„^> {- ^} < 27) 

Now letting r , o — > such that o" 2 /r — > 2D we obtain ($ = p = p(x, £)) 



| = £ GP ^ ^ (28) 

= e -"*' 2 Uj^J e "* /2 P-Pe"* /2 (^) ^"* /2 . (29) 

i.e. equation (2). Therfore, the dynamics of a particle on a fully extended 
polymer is governed by ordinary diffusion. 

Let us now return to a flexible polymer performing conformational changes. 
In this case the geometrical factor in the rate (21), f(x — y) follows a power- 
law [9]. Thence, we investigate the case of long-tailed transition rates. Let 



C^\x\' {1+ ^ \x\ > e > 



(30) 
otherwise 
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with = T(l + (j,) sin(7r/i/2)/7r and < /i < 2. Here e is the minimal jump 
length, which substitutes a fixed jump length a of our previous example. Note 
that we do not require f(x) in (30) to be normalized to unity. Any prefactor 
independent of e occurs in both terms in the master equation (19) and can 
be absorbed as a time constant. The particular choice of C M will become clear 
below. Inserting (30) into (21) gives 



dp(x, t) 
dt 



-/3*(x)/2 



dy 



\x-y\>e 



\x — y\ 1+ ^ 



-e^' 2 p{x,t) 



dy 



\x-y\>e 



e -/»(i/)/2 
\x — y\ 1+fl 



In the limit e — > (31) becomes (for notational ease we set s(x) 
q(x,t) = p(x,t)/s(x)) 



(31) 



= e -/3*(x)/2 



dp(x, t) 
dt 



= lim C u < 



s(x) 



dy 



q(y,t) 



\x-y\>e 



x — y\ 1+fl 



\x-y\>e ) 



lim C a < 



s(x) 



dy 



\x-y\>e 



q(y,t) - q(x,t) 
\x — y\ 1+tl 



-q(x,t) J dy 



\x-y\>e 



s(y) - s(x) 
\x — 



= CAs{x) 



dy 



q(y,t) - q(x,t) 
\x — y\ 1+tl 



(32) 



(33) 



(34) 



Note, that in the limit the resulting rates w(x\y) are not normalizable. There- 
fore, as mentioned earlier, the process is not a pure jump process in this limit. 
Implicitely, in the limit carried out above, not only the minimal step size e but 
also the typical waiting time r w (y) = (j dxw^y)) -1 at position y vanishes, 
because of the singular nature of the rate w(x\y) in this limit. However, the 
limit can be carried out since the resultant singularities cancel and eq.(34) can 
be interpreted consistently. 



The integrals appearing in (34) are symmetric fractional generalizations of the 
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ordinary Laplacian (see Appendix 1) 

At = r(u^ W 2) / /fa>-{W , (35) 

7r 7 |x — y\ ^ 

and equation (34) can be recast into a more concise form: 



^ = e-^' 2 At e^ 2 p - pe^ 2 At e~^ 2 (36) 
= ^p. (37) 

Again, the similarity to the underlying master-equation is obvious. Trivially, 
detailed balance is retained in the limit and the Boltzmann distribution is the 
stationary state of (36). Note that for ji = 2 we recover system (28), i.e. 



-G,2 — l-FP- 



(38) 



However, in the range < /i < 2 the operator £g,v does not portray a situation 
equivalent to any of the equations discussed above and does not correspond 
to any Langevin equation of the type 

dX = F(X)dt + dL^t) (39) 

in which the long-tailed influence is incorporated into the thermal fluctuations 
by a symmetric Levy stable process L^t). In [37] the difference between the 
fractional Fokker-Planck equation corresponding to Levy-stable additive noise 
and the topological generalization is investigated in detail. 



4 Comparison of systems in a double well potential 



We proceed by showing that the two generalizations of the FPE leading to 
Boltzmann statistics in equilibrium (one based on the fractal time, another 
on the fractal space approach) differ strongly in their predictions concerning 
relaxation in the same kind of potential. For our comparison we have chosen 
a classical Kramers situation, i.e. diffusion in a double-well potential. As a 
generic double well we chose 

= x 4 - 2x 2 + 1 (40) 

depicted in the upper left panel of figure 3. Initially a particle is placed at the 
potential minimum at x — — 1, p(x, 0) = 5(x + 1). The solutions p(x,t) are 
depicted at three different times (t = 0.03, 0.1, 1.0) in figure 3 for a value (3 = 
5.0 and an exponent [i — 1. The fractional diffusion equations were numerically 
solved by first mapping them onto the corresponding hermitian Hamiltonian 
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Figure 3. Evolution of densities p(x,t) in the double well potential <3?(x) (upper left) 
at three different times: t = 0.01 (upper right), 0.1 (lower left) and 1.0 (lower right) 
for (3 = 5.0 and fi = 1 with initial condition p(x,0) = 5(x + 1) in relation to the 
stationary state (thin solid line). The diffusion process equilibrates quickly within 
the left potential well and does not pass the barrier for short times. In contrast, the 
super diffusive processes transfer probability even for small times and remain peaked 
at x = — 1 for much longer. 

H and subsequent spectral decomposition TL = J2n ^nP n where the A n and V n 
are eigenvalues and projectors onto the corresponding eigenspaces. Solutions 
p(x, t) were then expressed in term of these quantities. 

Let us first note that the time-scales of the relaxation processes for the case 
of normal diffusion and for Levy-processes differ vastly. The typical times 
of the concentration equilibration (being the inverse of the largest nonzero 
eigenvalue) are: tq = 20.36 for the diffusion process (2), tq = 4.51 for the 
subordination case (/i = 1) and tq = 3.93 for the case of long jumps with the 
same value of [i — 1. A much greater value of Tc for the ordinary diffusion 
process is not surprising. In a short time it equilibrates within the left potential 
well. However, due to the continuity of its sample paths, the diffusion process 
accumulates probability in the right potential well only for long times. In 
comparison, the subordination process as well as the process governed by (37) 
can pass the potential maximum at x = even for short times, a direct 
consequence of the possibility of initiating a long jump from left to right. 
The possibility to transmit probability through the potential barrier is payed 
off by a slower equilibration within the left potential well in which p(x, t) for 
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Figure 4. The shape of pi{x,t) (left) and pji(x,t) (right) for subordinative and 
topological super diffusion at t = 0.1, (3 = 5.0 and exponent ji = 1. Note the difference 
in scale on the left and right, respectively. Despite the fact that these processes are 
identical if $ = 0, they respond quite differently if a potential is imposed. 

both, subordination and topological superdiffusion stay sharply peaked around 
x — — 1. This cusp-like shape in the left potential well is reminiscent of the 
classical Scher-Montroll picture of CTRW. Despite their qualitative similarities 
in the overall shape of p(x, t) subordination and topological superdiffusion are 
quite distinct on close inspection as is shown in figure 4. 



In order to avoid huge difference in scales, we plot in figure 4 the reduced 
densities in the left and in the right potential wells, i.e. the densities which 
are normalized to unity on the given sub-interval and for the given time: 



"^zdSfe <41) 

<42) 

Clearly, the shapes differ strongly on both sides of the potential. 

In order to reveal the differences in barrier penetration it is instructive to 
consider the time dependence of the total probability in the right potential 
well 

oo 

Q(t) = J dxp(x,t), (43) 
o 

and the rate r(t) defined by 

Asymptotically, as t — > oo the processes investigated here evolve according to 
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Figure 5. The time dependent rate r(t) for ordinary diffusion, subordination and 
topological super diffusion in the double well. Parameters are (3 = 5.0, fj, = 1. 



p{x, t) « p s (x) - (p s (x) - po(a:)) e Al * A x < 0. (45) 

where p s (x) and po( a; ) denote the stationary and initial states of the system, 
respectively. Therfore, r(t) — > — Ai as t — > oo, the rate approaches the inverse 
of the relaxation time constant of the process. In other words, the function 
r(t) measures how fast the probability transfer across the barrier approaches 
the quasi-equilibrium value. For the same parameters as in figures 3, 4 the 
rate r(t) is depicted in figure 5. The ordinary diffusion rate is initially zero 
and increases as soon as the process equilibrates within the left potential 
well. It then levels off to a relatively small constant value proportional to the 
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probability flux across the barrier. In contrast, for the superdiffusive systems 
the rate to transfer probability over the barrier is finite even for very small 
times. For the latter r(t) increases more slowly than in ordinary diffusion. 
However, the asymptotic rate is higher. Note that the subordination process 
penetrates the barrier at relatively high but almost constant rate for the entire 
time, whereas for topological superdiffusion the rate is initially smaller but 
increases for large times beyond the subordination rate. 



5 Conclusions 

We discussed different generalizations of a Fokker-Planck equation for the case 
of Levy processes in external force fields. We show that such generalization 
can proceed along different lines (assuming strong temperature or mobility 
fluctuations or allowing for long jumps), which lead to different forms of the 
generalized operator. We confined ourselves to situations in which the corre- 
sponding equations are thermodynamically sound and their solutions approach 
Boltzmann equilibrium. The different relaxation behaviors of Levy-processes 
were investigated and compared for the paradigmatic case of a double-well 
potential and deviations from normal diffusion were clarified. 
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A Fractional Differentiation 

The operators involved in fractional differentiation are generically non-local. 
Thus, depending on the particular choice of boundary values a variety of gen- 
eralizations exists, all of which suit a specific need in physical applications. 
Depending on the specific problem at hand different representations of frac- 
tional differential operators may be chosen. Here we concentrate on a fractional 
generalization, A M//2 , < fx < 2 of the ordinary Laplace operator A generi- 
cally encountered in superdiffusive systems. The reader is refered to [38] for 
additional information on fractional calculus in general. 
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A.l Definitions 



Incidentally, the most natural way to introduce fractional differential operators 
is based on the generalization of the n-fold iterated integral of a function f(x) 

X X 

= T^W / dyf{y){x-vT- 1 = — J dyf{y){x-yT~ l (A.l) 

a a 

to arbitrary fractional order a > 0, 



X 

ai:f(x)^-^- [ dyfiy^x-yy- 1 (A.2) 



r (a) 



a 



J:f(x) = ^-r I dyf(y)(y-xr-\ (A.3) 



x 



where we distinguish between the right- and left-handed fractional integrals 
a I x an d x Ia j respectively. A fractional derivative of order /i which is frequently 
employed in physical applications is given by the n th ordinary derivative of the 
fractional integral of order < n — fi < 1, where n = [/x] + 1. In other words, 



aD * f{x) = ~dx^ J: ~" (A - 4) 
d n 

* D « f{x) = d^ Jr "- (A - 5) 

Note that since the distinction between left- and right-handed fractional in- 
tegrations is passed on to fractional differentiation, fractional derivatives are 
not symmetric with respect to interchanging their subscripts. The definitions 
(A. 4, A. 5) imply that 

d m 

n m = -f— = (-l) m x D™ m = 0,1,2,... 

dx m 

That is, fractional differentiation is proportional to ordinary differentiation for 
integer values of the exponent. Q By sequential partial integration (A. 4) can 
be recast into 



1 To recover symmetry for integer values of the exponent m one frequently en- 
counters an alternative definition of the left-handed fractional derivative, namely 
x Daf(x) = (—l) n -^T X Ia~^ which includes the pre-factor (— l) n in the definition. 
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Therfore, the operations of fractional differentiation and integration generally 
do not commute. However, if we let a — > ±00 and require that lim^-too p k '(a) < 
00 (A. 4) and (A. 5) yield 



1 7 f {n) (v) 

,D£f(x) = —, r / dy . J ^ (A.7) 

xJ K ' r(n-//) J (x-yY +l - n v ' 



-00 
00 



1 7 f (n) (v) 

x D^f(x) = —^ / dy , J ;Z • A.8) 

In fractional diffusion equations the symmetric fractional generalization of 
the Laplace operator frequently appears. Its one dimensional variant can be 
defined in terms of a D£f(x) and x D£f(x) by 

A*ee-- ^--{-00^+,^} 0</i<2. (A.9) 

2 COS(7T / U/2) 

The particular case of requires some care. Although the definition (A.9) 
does not give a useful result for /x = 1, since both numerator and denominator 
vanish, the limit \i — > 1 can be interpreted consistently as will become clear 
below. Partial integration of (A. 7, A. 8) and insertion into (A.9) yields 



2r(-/i)cos(7r/i/2) J tf \x-y\ 1+ » 

— OO 

Since lim AA _,ir(— //) cos(7r/x/2) = — tt/2 the operator A2 is well defined, i.e. 



( Jf(y)-m \ 

it J * (x-y)< 



= r / dy UK "! T - (A.11) 



Clearly, A 5 ^ V = d/dx, the "square-root" of the Laplacian is symmetric, 
whereas the first derivative is not0The operators -oo-D£? xD^ and A 2 , possess 
a particularly simple form when the equation involved are Fourier transformed. 
If we define the Fourier transform of f(x) by 

T[f]{k) = f(k) = J dxe ikx f{x) (A.12) 

the transformed equations (A.7), (A.8) and (A. 10) read 



In the literature the symbol is frequently used instead of A 2 with can lead to 
some confusion if [i = 1. 
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F[. 00 D»f](k) = (-ik)y(k) (A.13) 
^D&f\(k) = (ikyf(k) (A.14) 
r[ASf](k) = -\h\»f{k). (A.15) 

Thus, the Fourier transform of fractional differentiation is equivalent to a mul- 
tiplication by Because of this, the investigation of equations containing 
operators such as A 2 is usually much easier when expressed in Fourier space. 
The simplicity of eqs. (A.13), (A.14) and (A.15) is often a reason to chose the 
latter as definitions of fractional derivatives. On the other hand, the corre- 
sponding integral representations (A. 7), (A. 8) and (A. 10) offer a more direct 
picture of the characteristic properties of the corresponding operator. For ex- 
ample, the Cauchy process obeys d t p(k,t) = — \k\p(k,t) in Fourier space, a 
simple equation indeed. By inverse Fourier transform this gives 



d t p(x,t) = (A*p)(x,t) (A.16) 
dy [w(x\y)p(y, t) - w(y\x)p(x, t)} (A.17) 



w(x\y)= 1 (A.18) 
tt(x — y) 2 

In contrast to the Fourier representation it is obvious from eqs. (A.17) and 
(A.18) that the fractional operator A2 is non-local, involves singularities and 
appearing in the example given above describes a jump process, a fact not 
readily available in Fourier representation. 

To complete the picture, let us mention that the form (A.15) of the fractional 
Laplace operator can be generalized to arbitrary dimensions n. Let /(x) be 
a function on M. n and /(k) the corresponding Fourier transform, then A2 
defined by 

^A*/](k) = -|k|"/(k) (A.19) 
yields upon Fourier inversion the integral representation 

^ 2/T (^) 7 , n My) - /(x)] 



Af/(x) = / d V V . I (A.20) 

y ' W 2 r |x-y|M+« 1 > 
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